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ABSTRACT: 



This is the third year-end report on the EHD contract. 

Work with the steam injector continued during the past year. 

The definition of the slip parameter for unsteady flow has been 

used in conjunction with velocity profiles measured with a hot wire 

anemometer to evaluate the potential role of EHD (or EGD) in anemometry. 

A direct measurement of particle size (using holographic interferometry) 

-4 

indicated that there are no particles of size 10 m or greater. A 
computer program describing the EHD flow has been written. 



FOREWORD 



This is a continuation of the EHD Research program started on 
September of 1968 at the Naval Postgraduate School, Monterey, California. 
This work was sponsored by the Naval Air Systems Command under the 
technical cognizance of Dr. H. R. Rosenwasser. 

One Master thesis was generated and it is included in the Appendix. 

A paper entitled "Unsteady-Flow Electrogasdynamics" is presently being 
submitted for publication. A Navy patent has been applied for the 
spectral anemometry principle. A progress paper was given at the 
Energy Conversion Colloquium held during 9-10 November at the Naval 
Postgraduate School. 

The help of LCDR Bohley and Mr. Patrick Hickey is gratefully 
acknowledged. The laser holographic measurements were made possible 
through the assistance and cooperation of Professor Collins of the 
Aeronautics Department, NPS . 
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I. INTRODUCTION 



During the past year progress was made basically along three fronts: 
First, the anemometry application was the subject of an LG thesis 
(included here as the Appendix). Also, a paper tying-in the work on 
anemometry and the EHD generator was written and is presently being 
submitted for publication. Next, an attempt was made at measuring 
particle size directly (i.e., by laser holography). Last, a computer 
program which solves the combined Laplace -Poisson equations inside 
the generator was written and checked-out. 

This report includes the final hot-wire anemometer frequency profile 
of the region being investigated. As pointed out in last year’s 
report"*", these measurements were made along the centerline in the wake 
of the cylinder at a location 5mm downstream of the exit plane of the 
injector nozzle. A review of the optimum frequency response concepts 
follows the hot wire data. Also included are a discussion of the 
laser -holograms and a brief summary of the computer program. 

II. HOT-WIRE ANEMOMETER READINGS 

The unit that gave the final results reported here is a Thermo 
Systems, Inc. model 1051-2 constant temperature hot wire, with a 
frequency response of 50 KHz. 

The hot wire readings were to be used as a standard of comparison 
for the EGD data but because the hot wire calibration is not appropriate 
for two phase flows, air was substituted for the steam in the cylinder'*'. 
The air flow was matched so as to give the same momentum as the steam. 
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The absence of droplets extended the life of the hot wires which was 
very marginal to begin with due to the intense vorticity of the wake 
of the cylinder. Another complication arose from the fact that the 
air flow through the steam nozzle was typically as a different temperature 
from the main air flowing around the cylinder. Several attempts were 
made at equalizing the two temperatures. The most reliable data 
seemed to indicate that the nozzle flow does not have a measurable 
effect on the main flow. This holds for distances of at least 5mm 
away from the nozzle exit plane. The resulting profile is shown in 
Figure 1. It represents the main air and the nozzle air (8 psig) 
but it is indistinguishable from the hot wire profile for the main 
air itself. There are two dominant frequencies apparent from this 
figure, namely, 1.2 KHz and 10 KHz. The first one represents twice 
the Strouhal frequency for the cylinder and is appropriately so 
because the measurements were taken along the centerline of the wake 
where the vortices , shed off the cylinder on alternating sides , add 
up. The Reynolds number based on cylinder diameter for this flow 
(55 m/sec) is 8x10^ so that the boundary layer on the cylinder is 
laminar. The second important frequency apparent from Figure 1 
seems to be associated with the nozzle that injects steam into the 
wake, see Figure 2. The 10 KHz peak appeared whether or not there 
was air issuing from the nozzle. The exit diameter of this nozzle is 
about 2mm which is some ten times smaller than the cylinder diameter 
and perhaps appropriately the frequency is nearly 10 times higher. 

It was observed further that 10 KHz is apparently a high enough 
frequency so that the dissipation of the vortices is high; at a station 
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10mm away from the nozzle exit plane, the 10 KHz peak disappears 
whereas the 1.2 KHz peak remains nearly as intense. So the fluid 
viscosity can act on these "lesser whorls" but the flow is still to 
be considered pseudo-turbulent since the vortices at the Strouhal 
frequency are basically two dimensional. 

There are no other significant contributions to the vorticity 
because the flow turbulence is carefully removed upstream and al 1 
measurements are made away from the wall boundary layers . There is 
one uncertainty that remains, however, and that involves whether or 
not the presence of large particles (i.e., two-phase flow) shifts 
the intense part of the profile towards the lower frequencies. The 
larger particles can absorb some of the eddy energy more or less 
continuously at this may cause a significant shift. This could be 
described with a different "effective" viscosity for the two-phase 
flow. The only way to resolve this dilema is to make measurements 
with either a laser or an ultrasonic anemometer. The laser device 
relies on the presence of large particles and measures their Doppler 
shift. The accuracy of these devices is considered below that of the 
hot wire at the present. 

III. OPTIMUM FREQUENCY RESPONSE 

It has been pointed out"*" that in the anemometer application of 
the EGD probe the charged particle size determines the frequency 
response of the device. This is clear insofar as the inertia of the 
particle governs how fast it can equilibrate with the eddies . 
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In last year's report, it was shown that the slip parameter (S) 



for unsteady flow could he written as 

max V(cu) + T v ^ 



( 1 ) 



where tu = frequency 

p = particle mobility 

E^ = breakdown field strength 

V(u>) = RMS intensity of the free stream velocity at the 
frequency tu 

t = equilibration time for the particles 

Since the mass of the particle only appears in t , the question arises 
as to what may be a desired criterion for the frequency response as a 
function of size which si consistent with a small slip parameter. This 
criterion maybe arbitrary, i.e., fixing iot to be some small fraction. 

Figure 3 shows a plot of a> vs radius for values of out equal to 

-2 -3 

10 and 10 . This choice fixes <u without regard to the value of 

_2 

the first term in Equation 1 which in itself may be higher than 10 ; 

however, a high slip due to a high mobility is undesirable for either 
the generator or the anemometer. Note further that if V(tu) is very 
small i.e., close to zero, then the slip may be very high but since 
there is no significant contribution by the flow at this frequency the 
slip value is unimportant. 

Since the slip parameter plays a significant role in the generator 
efficiency, it is necessary to determine the spectral distribution of 
the flow inside the generator and calculate the slip at each of the 
dominant frequencies. 
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IV. LASER HOLOGRAPHY 



In last year’s Final Report 1 , results from charge distribution 
measurements seem to indicate the existence of charged droplets of 

g 

size of about 100 p m. This is because about 3 x 10 electrons per 

pulse had been measured at a steam pressure of 13 psig. The resolution 

that we could obtain with the laser apparatus was itself about 

2 •? 

100 pm and the limiting factor ’ was the graininess of the hologram 
reproduced object. As can be seen from Figures 4 and 3 , no particles 
can be detected and we conclude that the particles are of a size much 
smaller than 100 pm. Note that the nozzle exit diameter is 0.l4 
inches or 2.5 mm. 

Other work in last year's Final Report seem to indicate that sizes 
were in the order of 10 pm and less. The direct measurement of these 
particles is a the limit of the technique even with a refined set-up. 

V. COMPUTER PROGRAM 
4 

As shown by Minardi the potential distribution within an EGD 
channel is described by a "combined Laplace-Poisson" equation. The 
two equations are required by the two-phase space charge nature of 
the EGD environment. Poisson's equation models the space-charge 
region, while Laplace's equation models the space-charge free region. 
Electrode geometry and space-charge effects on the potential distribution 
are of interest in the design of EGD devices. The solution of the 
governing equations will show these effects. However, a closed-form 
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solution is often very difficult to find, and any attempts at a solution 
will be complicated by the varying boundary conditions imposed by the 
electrodes. 

In order to have an analytical input to the design process, a 
numerical solution method (successive over -relaxation) was programmed 
for the IBM 360 computer system. Finite difference equations, derived 
from a Taylor-series expansion about the point in question, were used 
to build the numerical model. The program has the capability for 
solving axisymmetric and two-dimensional channel flows. The program 
output is designed to give both numerical and graphical representations 
of the solution. 

Of importance in the computer modeling of the EGD device is the 

5 

fluid dynamics involved. By means of a similar program, Minardi 
has shown that for particle sizes of interest, the fluid dynamics may 
be uncoupled from the electrodynamics. Thus it is possible to limit 
modeling of the dielectric gas flow to a constant rate of jet expansion 
and insure conservation of charge within the jet. 

Input variables to the program typically include those describing 
the electrode geometry, dielectric characteristics, charged particle 
characteristics, flow control, variables matrix dimensions and the 
output desired. The program simulates a steady state picture of the 
entire EGD situation by the superposition of two separate matrices. 

One matrix will carry all potential information, while the other matrix 
carries both charge distribution and some flow control information. 

The initial part of the program is devoted to setting up these 
matrices . 
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Once this has been accomplished, the program calls for an iterative 
calculation to determine the potential distribution, and a calculation 
for the |E | -field distribution. Remainder of the program is devoted 
to output of these results . 

The program was applied to the geometry shown in Figure 2. The 

dielectric gas used was steam issuing at 397 m/sec from the nozzle. 

This corresponds to a reservoir conditions of 8 psig and 2Uo°F. This 

is a slightly super-heated steam. Particle charging was considered to 

-7 2 

be corona limited, giving a particle mobility of 1 x 10~ m /v-sec. 

-12 

A dielectric constant of 8.9 x 10 coul/v-m is used. Attractor 
potential was 3 KV, and the collector potential was 50 volts below 
ground. Figure 6 is the plotted output of the potential distribution, 
showing lines of constant potential. Figure 7 is the plotted output 
of |E | - field distribution, showing lines of constant value of E. 

The potential has been normalized with respect to attractor or corona 
potential, and the E-field values with respect to the breakdown field 
of 3 x 10^ v/m. The geometry is represented by the cross-hatched 
points . 



VI. CONCLUSIONS AND RECOMMENDATIONS 

It is evident from the results shown in the Appendix that we have 
been able to achieve "acceptable" responses from the EGD probe up 
to 1000 Hz. This acceptability is based on the correctness of the 
hot-wire data used for comparison. From Figure 3 S this implies that 
particle sizes are centered about 1 pm or higher radii. From the 
holographic measurements, we know the particles are not much above 
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10 (i,m. Our efforts at changing the superheat of the steam in order 
to affect particle size have indicated that this is not the only- 
relevant parameter; apparently what happens to the charged particles 
outisde of the nozzle is equally important. There, of course, the 
behavior is somewhat random and a statistical sample is needed. 
Lowering the pressure from 13 psig to 8 psig is a step in the right 
direction but apparently something happens that gives multiple peaks 
in the distribution of particles . These results point out two things 
First, better humidity control and some dilution of the steam inside 
the cylinder (through, say, a bleed of nitrogen) may be needed to 
achieve a reduction of particle size. Second, the present response 
is sufficient to measure some interesting turbulent flows which need 
not be as complicated as the wake of the cylinder (for example, 
grid-generated turbulence). 

The computer study will be directed toward the design of guard- 
electrodes inside the generator channel which will improve the output 
of the device. The program is further useful in pointing out regions 
of high electric field. This is quite significant for any EHD power 
generator . 

The definition of the coupling of the charged particles with an 
unsteady flow field remains as an important goal of this project. 

Most EHD flows are unsteady to some extent and a low slip parameter 
for steady flow may not be good enough if significant velocity 
components exist at non-zero frequencies. 
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FIGURE 2 EGD GENERATOR SCHEMATIC 
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FIGURE 3 OPTIMUM FREQUENCY RESPONSE 
AS A FUNCTION OF THE RADIUS 
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FIGURE 4 LASER HOLOGRAM OF VIRTUAL IMAGE 




FIGURE 5 ENLARGED VIEW OF TEST SECTION 
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FIGURE C TWO DIMENSIONAL OUTPUT OF POTENTIAL DISTRIBUTION 
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FIGURE 7 TWO-DIMENSIONAL OUTPUT DISTRIBUTION OF 
THE ELECTRIC FIELD 
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ABSTRACT 



This study involves the measurement of RMS velocities 
covering the frequency spectrum 25 to 20,000 hertz using 
an electrogasdynamic (EGD) probe. Previous work proved 
the feasibility of this investigation. The flow of air 
behind a circular cylinder was measured. The freestream 
velocity was 185 ft/sec. Charged water droplets entering 
the wake through a jet in the cylinder at sonic velocity 
acted as the media for conveying information to the probe. 
Hot-wire anemometer readings made in the same spectral 
range were used as the standard of comparison. However, 
these need further refinement. 

Small droplet formation was optimized by determina- 
tion and control of the relative humidity of the freestream, 
steam temperature and steam pressure. Optimum droplet 
charge was then determined and controlled with corona 
voltage. A 10 to 15 decibel increase in high frequency 
response and stable performance resulted. 
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I. INTRODUCTION 



Previous work in Electrogasdynamics has indicated a 
definite correlation between EGD probe readings of col- 
lector current and the mean and rms velocities measured 
by a hot-wire anemometer in the turbulent wake of a cir- 
cular cylinder [6]. This thesis is an attempt to improve 
the spectral correspondence of EGD probe measurements to 
hot-wire results. The probe is a stainless steel rod with 
a sharpened tip which is placed in a moving stream of air 
where it collects charged particles. For charged particles 
of proper size and charge, sufficient coupling may result 
in the turbulent flow to allow measurement of the spectrum 
of turbulence intensity. 

The advantages of the EGD probe as a diagnostic tool 
over the hot-wire anemometer are manifold. The probe is 
very rugged, allowing probing of intense turbulent fields. 
It is very simple and inexpensive and, therefore, can be 
easily manufactured. The probe tip area is small and sym- 
metrical, allowing better spatial resolution than the hot- 
wire. Since the very small charged particles can react to 
high frequency oscillations, the frequency response of the 
EGD probe may be higher than that of the hot-wire. In ad- 
dition, signals from the EGD probe, which are electrical 
in nature, need little or no further amplification [7]. 
Effects corresponding to the thermal inertia of hot-wires 
are not present. There is no requirement to feed the 
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output through a linearizer since the EGD output is already 
linear, i.e., a direct measurement. Also eliminated is the 
effect that free convection currents have on the hot-wire 
zero setting. In the case of two-phase fluids, the EGD 
probe needs no separate calibration and will withstand 
the impact of charged droplets or particles whereas the 
fragile hot-wire will readily break under these conditions. 
The main disadvantages of the EGD mechanism are the require- 
ments for a source of charged particles of optimum size and 
charge and for a high voltage supply (along with the neces- 
sary safety precautions) . The optima appear to be rather 
critical . 

The specific objectives of this thesis are: 1) im- 

provement of the high frequency response of the EGD probe 
[5] ; 2) the elimination or explanation of inconsistencies 
between the EGD results and the hot-wire results; and 3) 
the determination of the effect that various parameters, 
i.e., humidity, steam pressure and temperature, the ratio 
of collector current to total current, etc, have on the re- 
sults. A lesser objective was finding a means of facili- 
tating the reduction of data for comparison to hot-wire 
data . 

This work complements part of a comprehensive, long- 
term project being conducted by Assistant Professor Oscar 
Biblarz, Department of Aeronautics, Naval Postgraduate 
School. Existing EGD equipment installed in Building 230 
at the Naval Postgraduate School was utilized. 
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II. THEORY 



In 1959 Bennett [2] pointed out that by carrying elec- 
trical charge in the volume of a stream of air instead of 
on the surface of a belt, a gaseous Van de Graff generator 
could be built. Since a greater volume of charge would be 
carried at higher speeds, much higher current should re- 
sult. This marked the beginning of the EGD generator. 

The basic EGD mechanism involves particles which are 
introduced into an air stream or other suitable gas at the 
entrance of the test section (Figure 1) . For the source 
of charges an aerosol flow in conjunction with an ionizer 
are used. Saturated steam is injected through a corona 
discharge into the "turbulent" wake of a cylinder. The 
corona discharge is produced by applying a high voltage 
between a sharp needle and an attractor ring which is em- 
bedded in the wall of the steam nozzle located on the 
downstream side of a right circular cylinder. Electrons 
from the negative needle drift toward the attractor and 
attach themselves to the then-growing water droplets form- 
ing negative charged particles. Some of these charged 
particles are in turn dragged or convected by the air flow 
down the conversion channel where their electrons are 
transferred to a metallic collector probe. Since the col- 
lector is also negatively charged, work is accomplished 
and some of the kinetic energy of the gas has been changed 
directly into electrical energy. These charges flow 
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through an outside circuit and return to the needle. If 
the energy drained from the gas is minimal, then the flow 
has not been disturbed by the conversion process. 

One of the important parameters here is particle size. 
If the particle is too small, i.e., very mobile, the drift 
velocity of the particles becomes large and current loss 
results as exemplified in the current density equation for 
EGD, 

j = p e (U. - yE) (1) 

where yE is the drift velocity. By rearranging the equa- 
tion in its scalar, one -dimens ional form one obtains 

j = p e o-f>". t 2 ) 

00 



and slip can be defined as 



s 




(3) 



For a certain range of radii, as the charged particle size 
decreases, the mobility constant increases (Figure 2) and 
hence slip increases [3,9]. 

Mobility is a function of mass and charge of the 
particle. Since the ability of a charged particle to fol- 
low the turbulent eddies in the cylinder wake is a function 
of its mass, an increase in size results in lower frequency 
response of the charged particle. This can be seen by con- 
sidering the effects of inertia in a moving particle. A 
particle with less mass can accelerate, decelerate, and 
change direction much faster than a particle of greater 
mass exposed to the same forces. 
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One of the problems in previous investigation [6] 
has been the generation of steam droplets of correct size 
for high-frequency coupling and correct charge for low 

mobility. Injection of particles ranging between 10 ^ 

- 7 

to 10 meters charged in a corona field results in a 
maximum drift velocity of one to ten meters per second 
but allows particles to fluctuate with frequencies up to 
50 kHz [4]. Previous work showed a high frequency drop 
for the EGD probe in comparison with hot-wire anemometer 
measurements [6]. An analysis indicated that the particle 
radius was on the order of five microns [4] . Therefore it 
seems that the size of the steam droplets can be reduced 
by decreasing the condensation which should, in turn, re- 
sult in an increase in EGD probe bandwidth. Previous work 
showed that "condensation" was decreased by lowering the 
steam pressure to below 8.3 psig [4]. Wetness could be 
decreased by superheating the steam as it travels from 
the steam generator to the cylinder exit nozzle. In ad- 
dition, Reference 13 (unpublished) and some of our pre- 
liminary work indicated that certain ranges of relative 
humidity of air in the channel could enhance particle size 
after injection into the air flow. Thus, high-frequency 
response could also be enhanced by controlling relative 
humidity . 

One last area concerns the standard of comparison. 

The hot-wire anemometer measures the turbulence in a 
single-phase flow while the EDG probe measures a two-phase 
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flow. The frequency spectrum of these two 
fer since the droplets most likely enhance 
quency spectrum. 



flows may dif- 
the low fre- 



15 



Ill . EQUIPMENT 



The heart of the apparatus is the small plexiglass 
wind tunnel test section. The inside dimensions are two 
by four inches in cross section and ten inches in length. 
The test section contains a steam injector unit which is 
fed by a heated steam line from a steam generator [10]. 

The aerosol injector nozzle consists of a stainless steel 
ring embedded in the cylinder wall which allows the in- 
jection of condensing steam into the turbulent wake of 
the cylinder (Figure 3) . 

The steel ring also serves as the corona ring posi- 
tive electrode of the aerosol ionizer. The negative 
electrode consists of a steel needle positioned on the 
corona ring centerline but slightly recessed. A Sorensen 
High-Voltage DC Power Supply provides the high potential 
between the ring and needle necessary to negatively charge 
the steam particles as they pass through the throat of 
the nozzle (at sonic velocity) into the turbulent cylin- 
der wake. The potential difference is monitored by a 
Sensitive Research electrostatic voltmeter. A few milli- 
meters downstream from the cylinder nozzle a third elec- 
trode, consisting of a stainless steel probe with a rounded 
tip (Figure 4) , collects the current and thereby detects 
the frequency spectrum of the turbulent wake via the infor- 
mation carried by the charged particles. 

The air is supplied by a Carrier three-stage centri- 
fugal compressor with a maximum flow rate of 4000 cubic 
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feet per minute and a maximum pressure ratio of two. This 
air is fed through a cooling bank into the main plenum 
which is attached to the plexiglass test section. The valve 
on the cooling bank which controls the amount of air being 
cooled also enables control of the temperature and hence 
the relative humidity of the air (Figure 5) . The air passes 
through the test section and around the cylinder, then 
mixes with the steam jet. The particles are convected to 
the collector probe, where they are neutralized, and the 
water-air mixture flows out of the test section into the 
atmosphere . 

The current from the collector probe is fed into a 
black box which houses the circuitry for the probe (Figure 
6) . The three International Rectifier diodes prevent re- 
verse current flow [11]. A Simpson microammeter measures 
the collector current. The capacitor and the neon bulb 
safety switch protect the frequency analyzer. The capaci- 
tance is selected large enough to ensure a phase shift of 
less than one degree in the frequency spectrum from 25 
hertz to 20 kilohertz and to accomodate the high voltages 
which are present (Appendix A) . 

The pre -amplifier (Figure 7) is a self-powered unit 
consisting of a UA741C operational amplifier, resistors, 
and a battery. It provides a gain of ten over a frequency 
range of 100 kilohertz. A shunt resistance of 39,000 ohms 
is placed across the input in order to avoid saturating 
the operational amplifier in the preamplifier, i.e., one 
third of the signal is bled off by this shunt. 
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The frequency analyzer is a General Radio Company Type 
1921 Real Time Analyzer (Figure 8) . This consists of a 
1925 Multifilter and a 1926 Multichannel RMS Detector. The 
signal (Figure 9) from the preamplifier is put into the 
multifilter, separated according to frequency, put into the 
RMS detector where the data are converted to binary digital 
form and the rms level computed by true linear integration. 
The computed data are then displayed digitally on the rms 
detector nixie-tubes and the analog form of the data is 
projected graphically on a Hewett-Packard X-Y display. 
Either a Tektronix Type RM 15 or a Type 551 Dual Beam 
oscilloscope is used to monitor the signal input to the 
frequency analyzer. 

The hot-wire data of this experiment were collected 
with a Thermo Systems, Incorporated, Model 1050 and 1051-D, 
constant temperature hot-wire anemometer (Figure 10) . The 
frequency response of the TSI equipment was superior to 
the frequency response of two other constant temperature 
hot-wire anemometers which were also tried, i.e., a Flow 
Corporation Constant Temperature Anemometer, Model 900, 
and a Security Associates, Model 100, hot-wire anemometer. 
In addition, the TSI unit had a built-in signal generator 
that enabled simplified procedures for frequency calibra- 
tion. A Kronhite Model 3202R bandpass filter was used in 
conjunction with the TSI unit to filter out the noise from 
outside the test frequency spectrum, which was 25 to 20,000 
hertz . 
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Other equipment included a psychrometer for measuring 
test section air humidity, a Simpson microammeter for 
measuring needle current and a voltage regulator for 
regulating line voltage. 
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IV. PROCEDURES 



The experimental procedure developed from an initial 
familiarization with the operation of the apparatus and a 
desire to start with a verification of previous work. This 
produced results similar to those reported in Reference 6. 

The succeeding runs involved a change from a 0.5 
microfarad, 2000 Vdc capacitor to a 10 microfarad, 2000 
Vdc capacitor in order to eliminate a phase shift in the 
portion of the spectrum below 250 hertz (Appendix A) . A 
closer correspondence in the lower half of the frequency 
spectrum now existed, between the EGD turbulence intensity 
(Tg) and the hot-wire turbulence intensity (T, ) , than 
measured in previous work [6], 

Asbestos insulation was wrapped around the bare por- 
tion of copper tubing immediately below the cylinder. This 
prevented early cooling, and subsequent condensation, from 
forming particles of too large a mass, which translates to 
a great inertia and a low mobility. During the test runs 
to check these changes, attention was given primarily to 
tape heat (this tape heats the line between the steam 
generator and the cylinder in the air stream) , steam pres- 
sure, and applied voltage. The temperature range was be- 
tween 230°F and 288°F, the pressure between 7.9 to 14.1 
psig, and the voltage between 215 to 3.1 kilovolts. 

Since much of the data taken on different days showed 
little correlation even when the above parameters remained 
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unchanged, it was necessary to investigate other parameters. 
Closer observation of the relative humidity, use of a more 
accurate .device for measurement of humidity, and measurement 
of humidity in the test section during operation revealed 
that the best results occurred on days when the relative 
humidity in the channel was less than 46%. 

For operation at 25% relative humidity, test section 
temperature ranged between 95°F to 125°F. Difficulty was 
encountered in getting values of collector current and 
total current greater than zero when the relative humidity 
was 25%. Even at 30%, Ig and I T tended to be erratic, the 
microammeters registering zero approximately half of the 
time. Good results occurred at 32% and this value was 
chosen for subsequent testing. 

The next portion of the experiment was devoted to a 
thorough testing of the equipment after it was discovered 
that the internal calibration voltage was occasionally er- 
ratic. The General Radio equipment tested in accordance 
with manual specifications. Noise in the output at the 
low frequencies (most prevalent around 60 hertz) was re- 
duced by grounding all equipment used in the experiment. 

A line voltage regulator was installed to decrease fluc- 
tuations in line voltage which might affect the equipment 
and results. 

A method of comparing existing hot-wire results with 
EGD measurements was devised using the multifilter unit. 

By subtracting out the values (in decibels) of the turbulence 



21 



intensity as measured by a hot-wire throughout the frequency 
spectrum on each of the multifilter channels, any deviation 
from a straight line on the X-Y display gave an immediate 
indication of the portion of the spectrum with which the 
EGD measurements did not favorably compare and the amount 
by which the intensity differed from the hot-wire data. 
Because it appeared that some values of the multifilter 
settings would subtract out a greater value of turbulence 
intensity than current EGD measurements were producing, 
the preamplifier was incorporated. 

The hot-wire anemometer portion of this experiment 
was conducted to obtain measurements of turbulence in- 
tensity, T^, at a nozzle pressure of eight psig and to 
verify previous data taken at 13 psig. In order to simu- 
late the jet effects but avoid the two-phase flow in which 
the hot-wire calibration is not useful and also to avoid 
breaking the hot-wires by droplet impact, the steam unit 
was disconnected and air at the same pressure was injected 
into the cylinder and out of the nozzle at 8 and 13 psig 
[10]. A limited amount of testing (which occurred before 
the failure of a field-effect transistor) with the Securitv 
Associates Hot-Wire Anemometer produced results which agreed 
with those of reference [6]. 

A Thermo Systems, Incorporated (TSI) constant temper- 
ature anemometer provided a much greater intensity at the 
high frequency end of the spectrum than had been previously 
measured (Figures 11 and 12) , especially when nozzle air 
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was on. Figure 11 was the result obtained using a DISA 
55F31 probe. All succeeding runs were made with probes 
used in previous work which are locally manufactured. 
Subsequent runs with a bandpass filter did not improve 
the results. It was later discovered, however, that a 
signal greater than 0.5 volts was possibly injected into 
the General Radio 1560-P40 preamplifier, damaging the 
field-effect transistor. Even though hot-wire measurements 
were made bypassing this preamplifier, it did obtain power 
from the frequency analyzer and was left connected to the 
analyzer. It is not known whether this could have caused 
a stray signal to be injected into the analyzer causing 
this unexpected, high frequency "tail." Another measure- 
ment with the TSI unit was attempted after the above ex- 
perience and culminated in failure to the TSI unit before 
any results could be obtained. 

The final portion of this experiment produced the 
EGD results reported herein. A smaller stainless steel 
probe (0.062 inch diameter) replaced the 0.125 inch steel 
probe used initially [6], It was felt that a smaller probe 
would sample a more localized area and not tend to average 
the turbulent eddies present in a larger volume of space. 

A new preamplifier (Figure 7) replaced the General Radio 
1560-40. A shunt resistance was incorporated to avoid 
over-driving the operational amplifier. Channel humidity 
was held at 32%, as measured at the exit of the test sec- 
tion, eight centimeters aft of the cylinder nozzle. Tape 
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temperatures were varied between 230°F and 295°F. Generator 
collector current was varied between zero and 10.2 microam- 
peres and applied voltage between zero and 4.43 kilovolts. 
Collector currents less than 0.2 microamperes were not mea- 
sured. 

The frequency analyzer converted the AC electrical 
input from AC voltage to the appropriate value in decibels. 
The measure of turbulence intensity, 7, at each frequency 
sampled, was expressed in dB and read directly from the 
X-Y display or the nixie tubes (Appendix B) . The pro- 
cedure developed for comparing the EGD data with the hot- 
wire results is as follows: 7^ represents the value of 

turbulence intensity measured with the EGD probe as given 
by the frequency analyzer and expressed in decibels. 7^ 
is likewise the value of turbulence intensity measured with 
the hot-wire anemometer. If the difference, A7 = T, - 7„ , 
is the same (± 1 dB for error) for each frequency, then the 
difference results solely from the difference in amplifica- 
tion of the two signals and the two signals are considered 
to be identical. In this context 7 = log-^ T^(f) (Appendix 
C). 

The procedure developed for determining the high fre- 
quency response involved a comparison of the difference in 
intensity between the most turbulent and the least turbulent 
frequencies. 7 is the maximum turbulence intensity, 
usually found at either the Strouhal or twice the Strouhal 
frequency. 7 . is the minimum turbulence intensity, in 
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decibels, normally found 'at the highest frequency. AT ucn = 

nrR 

T - T . is a measure of the difference in intensity 
of the turbulence between the two frequencies. If this 
difference is small, then a good relative high frequency 
response is indicated. ^-^hFR a l wa y s have a value 

greater than zero since, theoretically, turbulent vertices 
at high frequencies dissociate rapidly to lower frequency 
vertices [7]. The typical intensity versus frequency 
curve for a free jet exhibits a gradual decrease in in- 
tensity as the frequency increases [7] . The best (smallest) 
value of AT ucn obtained for both EGD and hot-wire data was 
12 dB. 



25 



V. RESULTS AND DISCUSSION 



Good results were consistently obtained when little 
or no tape heat was applied and when the collector current 
was very carefully controlled (Table I). Figure 13 shows 
a corresponding drop in high frequency response as temper- 
ature is increased. In order to examine and compare data, 
a plot of pressure versus density for various values of 
temperature (Figure 14) was constructed from steam tables 
[1,8]. This graph was used to convert values of pressure 
and temperature into corresponding values of density. A 
check of the steam pressure and temperature values imposed 
during the experiment, on this plot, reveals that best re- 
sults were obtained when points lay in the saturated or 
supersaturated regions. 

Two methods of comparing EGD and hot-wire data were 
used, both using the data shown in Figure 11 as the stan- 
dard of comparison. Figures 15 through 17 show typical 
results at 12 to 13 psig with varying temperatures and 
voltages. With the tape heat turned off, the temperature 
of the steam at this pressure is 240°F. Saturation tem- 
peratures of the steam for these pressures are 244°F and 
246°F respectively. Figures 18 through 20 show typical 
results at 8 psig at various temperatures and voltages. 

At this pressure the saturation temperature is 235°F and 
the steam temperature with the tape heat off is 230°F. 

These results are compared with the best hot-wire anemometer 
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data available. These data were chosen as the standard of 
comparison since these were the first data collected with 
the TSI unit and are expected to be free from any effects 
due to subsequent malfunctioning electronics. Since all 
data are already converted to decibels, to ratio (see Ap- 
pendix C) the desired points one simply takes the dif- 
ference of the logarithm of the plotted points at each 
frequency, measured from some arbitrary reference line, to 
the points to be matched. Twice the Strouhal frequency 
was chosen as a reference point since it is the most dom- 
inant point on the EGD and hot-wire data collected and 
also since the frequencies just beyond this point are most 
responsive to changes in parameters. 

The second method of comparing EGD data with hot- 
wire data was by inserting the hot-wire data into the multi- 
filter as explained earlier. Figures 21 and 22 show a 
comparison of typical results displayed during testing. 

Figure 23 shows the effect of the density of the 
steam on the current efficiency (I^/I^) for the two-phase 
fluid. In addition to the sensitivity of the current 
efficiency to a small change in voltage, the effect of an 
increase in density is to require a higher voltage in 
order to obtain the same current efficiency. This shift- 
ing of the curves to the right is similar to that reported 
by Whitby [12] for the effect of pressure (temperature was 
held constant) in a single-phase fluid (air). 
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Comparisons of Table I with Figure 13 and Figures 15 
through 22 reveal that the best high frequency response, 
i.e., when is least, generally occurs when the steam 

temperature is at or below saturation temperature for a 
given pressure. Figure 24 is a plot of pressure versus 
high frequency response for three ranges of temperature 
to saturation temperature ratios. The relative humidity 
and the collector probe distance remained constant. The 
corona voltage was held constant at three kilovolts with 
exceptions shown beside the data points. The best high 
frequency response (lowest ATjjp^) occurs at the super- 
saturation and saturation temperatures with the best re- 
sults in the eight psig pressure range. 

Figure 25 shows that a definite optimum voltage exists 
for operation at good high-frequency response. All other 
variables are held constant but a slight decrease in the' 
relative humidity from 34% to 32% and a slight increase 
in pressure from 8.1 to 8.2 psig during the experiment 
may have had an effect on the optimum collector current 
and on the current efficiency. The accompanying changes 
in collector current and current efficiency are shown. 

After a "triggering" voltage is attained, collector cur- 
rent increases linearly until finally a plateau is reached. 
The graph shows that the optimum current efficiency does 
not occur at the optimum voltage, although some previous 
observations indicated that good results occurred at a 
current efficiency of 0.2 to 0.5 and at a higher collector 
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current. Figure 26 shows the effect of current efficiency 
for three different temperatures with other parameters 
held constant. 
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VI. CONCLUSIONS 



It is concluded that an e lectrogasdynamic probe can 
indeed by successfully used to measure intense turbulent 
fields. Results comparable to or better than some of the 
existing hot-wire anemonmeters have been obtained. 

The following specific conclusions are drawn: 

1. A controlled relative humidity of 321 resulted in 
a significant improvement in the measurement of turbulence 
intensity in the high frequency range of the spectrum. 

2. Little or no tape heat gave optimum results (Fig- 
ure 22). This parameter is a highly critical one. A 
change of 1°F can, at times, make a difference of measuring 
a highly intense turbulence or not detecting any turbulence 
at 20 kHz. Sometimes saturation temperature provides op- 
timum results. At other times no tape heat (5°F below 
saturation temperature) provides optimum results. Even in 
this last case, an increase in tape heat of one degree re- 
sults in a marked steadiness of total current and collector 
current. An increase of an additional 1°F would steady 
the current even more but, for some runs, would also result 
in a significant high frequency detection loss. This fact 
is contrary to expectation, (earlier investigation indicated 
that the steam droplets seemed too large already) . 

3. A ratio of collector current to total current 
greater than 0.20 seemed to provide best results (Figure 26). 

4. A marked improvement could sometimes be seen dur- 
ing transient conditions, especially after a change in voltage. 
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The time required to reach steady state conditions is 
generally 30 seconds to one minute. However, in a few 
instances, transients were noted even after five minutes, 
indicating that steady state conditions are important. 

5. An optimum applied voltage range exists, below 
which no collector current is produced and above which 

no high frequency turbulence can be measured. This latter 
phenonemon may be due to an increasing charge built up 
per particle, resulting in an increased drift velocity 
back upstream to the attractor ring (Figure 25) . 

6. The effect of an increase in steam pressure is 
to shift the current efficiency curves to the right 
(Figure 23). This is comparable to conclusions drawn by 
Whitby [10] for single-phase flow. Saturation effects of 
the steam were lost in converting pressure and temperature 
to density, however, and are not accounted for in this fig- 
ure. These also seem to play an important role. 

7. Eight psig seems better than thirteen psig. This 
is in agreement with our previous arguments for smaller 
mass particles. 

8. Total collected current seemed to have an effect. 
For best results a collector current of one-half microam- 
pere or greater was required (Table I) . 
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VII. RECOMMENDATIONS 



Several areas need to be further explored in order 
to refine the electrogasdynamic probe and enhance future 
applications of this device. The areas are the follow- 
ing: 

1. Measurement of the turbulence intensity with 
state-of-the-art hot-wire anemometer equipment to deter- 
mine an accurate and consistent turbulence intensity fre- 
quency spectrum. 

2. An investigation into the effect of higher ap- 
plied voltages on the production of charged particles and 
a better high voltage supply, i.e., regulated to avoid 

± 60 volt fluctuations at 3.0 kilovolts (RMS ripple factor 
is 2.01). 

3. Better instrumentation for the control of steam 
pressure and temperature, i.e., more accurate monitoring 
of these data. Temperature control should be within at 
least ± 0.5°F, preferably ± 0.2°F and pressure control to 
at least ± 0.05 psi. This would allow accurate deter- 
mination of steady state conditions and insure that the 
critical parameters are indeed constant. Also, continuous 
monitoring of air humidity upstream of the cylinder is 
desirable to eliminate interruption of the experiment for 
humidity measurement and to preclude any unknown changes 
from occurring between measurements. 
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4. A study to determine whether there is any signi- 
ficant difference between the turbulent field of a single - 
phase fluid and that of the fluid when a second phase is 
introduced. It is possible that the two-phase flow may 
shift the spectrum towards lower frequencies by virtue 

of the large, massive particles which are present. 

5. Elimination of the temperature difference be- 
tween the free stream and the nozzle air. In addition to 
the lack of two-phase flow, using the hot-wire as a ref- 
erence introduces another discrepancy - the presence of 
two temperatures in the fluid. Since the hot-wire mea- 
sures temperature, any sensing of intermittent tempera- 
tures due to two different temperatures rather than to a 
difference in turbulence intensity could lead to a signi- 
ficant distortion of the frequency spectrum curve. 

6. A theoretical study of the nature and frequency 
range of the turbulence in the wake of the cylinder and 
nozzle to determine if, indeed, the hot-wire anemometer 
is measuring all effects present. 

7. Further investigation of the effect of small 
changes in temperature and pressure on the formation of 
droplets and the production of collector current at various 
frequency ranges in the turbulence spectrum. An accurate 
determination of the size of the droplets as a function 

of the various parameters, especially temperature and 
ppressure, at a given humidity, will greatly aid in deter- 
minimg the proper environment for the future use of the 
EGD probe . 
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TABLE I. Summary of Best EGD Results (L 



STEAM 

(FOR IONIZATION) 
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FIGURE I. EGD GENERATOR SCHEMATIC 
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Porticle Diometer {jj, m ) 

FIGURE 2. MOBILITY OF WATER DROPLETS IN AIR AT I* ATM. AND 20°C 
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FIGURE 4. 

COLLECTOR PROBE 
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AIR AND STEAM FLOW SCHEMATIC. 
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FIGURE 6. 

ELECTRICAL SCHEMATIC OF COLLECTOR 
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FIGURE 7. 

SCHEMATIC OF PREAMPLIFIER 



OUTPUT 
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X-Y DISPLAY, RMS DETECTOR AND MULTIFILTER 




MULTI FILTER UNIT 



FIGURE 8 

FREQUENCY ANALYZER 
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FIGURE 9. 

BLOCK DIAGRAM OF SIGNAL FLOW 








TSI HOT-WIRE ANEMOMETER AND FILTER 




Apparatus Ready for Testing 



FIGURE 10 

CONFIGURATION FOR HOT-WIRE TESTING 
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MAIN AIR ONLY (185 fps) 




MAIN AIR (185 fps) NOZZLE AIR (8 psig) 

FIGURE 1 1. X-Y DISPLAY FROM HOT _ WIRE ANEMOMETER 

L = 10 mm ( (£_) 

43 








MAIN AIR ONLY (185 fps) 




MAIN AIR (185 fps) NOZZLE AIR ( 8 psig) 

FIGURE 12. X-Y DISPLAY FROM HOT-WIRE ANEMOMETER 

L = 5 mm (Q_) 
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I g = 0.9/xA, MAIN AIR (185 fps), STEAM ( 8.0 psig , 231 °F ) 



Reference 



- 70 dB 

— 55 dB 
_ 40 dB 



FIGURE 13. X-Y DISPLAY FROM EGD PROBE 

L = 5 mm ( (£_) 
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FIGURE 14. PRESSURE -DENSITY-TEMPERATURE PLOT 
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15. COMPARISON OF EGD AND HOT-WIRE DATA 
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19. COMPARISON OF E6D AND HOT-WIRE DATA 
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FIGURE 20 . COMPARISON OF EGD AND HOT-WIRE DATA 
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AS ABOVE BUT WITH MULTIFILTER CORRECTIONS 



FIGURE 21 

COMPARISON OF EGD DISPLAYS 
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I = 0.6/iA, MAIN AIR (185 fps) , STEAM (8.0 psig, 230 °F) 




AS ABOVE BUT WITH 



MULTIFILTER CORRECTION 



FIGURE 22 

COMPARISON OF EGD DISPLAYS 
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FIGURE 24. OPTIMUM TEMPERATURE AND PRESSURE 
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FIGURE 25. OPTIMUM VOLTAGE AND CURRENT 
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CURRENT EFFICIENCY 



APPENDIX A: PHASE SHIFT CALCULATIONS 



The phase shift is calculated below in the order in 
which equipment changes were made. The relation governing 
phase shift is: 



2TrfCR 

1. C = 0.5 microfarads, R = 100,000 ohms. 
At f = 25 Hertz : 

-1 1 

<j> = tan r?r r 

2 (25) (0 .5 x 10 °) (10 s ) 



= 7.3° 

At f = 20,000 Hertz : 

<j> = tan" 1 (0.00016) 

= 0 ° 

2. C = 10 microfarads, R = 100,000 ohms. 
At f = 25 Hertz: 



= 0 ° 

3. C = 10 microfarads, R = 825,000 ohms. 
At f = 25 Hertz: 

4 > = 0 ° 





R 



1 

toCR 

1 



tan' 1 (0.127) 



(j) = tan 1 (0.0064) 
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APPENDIX B: DETERMINATION OF V, 



RMS 



The signals measured in this investigation are ex- 
pressed in decibels for ease in manipulation. The decibel 
is simply a ratio of some power level, W, to a reference 
power level, W , and defined as 



A zero dB gain implies that V = V Q . If V = 10V o then 

V DMC = 20 dB. Likewise, if V = 100,000V then V m . c = 100 
RMS ’ ’ o RMS 

decibels . 

The signal into the RMS detector is integrated over 
a four-second period during which each of thirty discrete 
frequencies are sampled simultaneously 1024 times in one- 
third octave bands. The detector measures values of vol- 
tage with a one -millivolt reference level. This is 
converted to decibel form using the equation 




o 




10 log 10 (V/V o ) 2 



20 log 10 (V/V o ) dB. 



N 




2 




V 



o 



dB 



i=l L 



where 



V(t^) = measured value 
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reference level 



V = 
o 

N = number of samples (binary number) . 

The computed data are then displayed digitally on the nixie 
tubes and are projected on the X-Y display in analog form. 

A scale-factor adder adjusts the output-level read- 
ings. Normally a one-volt RMS input produces an output- 
level reading of 60 dB . This reading is adjustable in one 
dB steps to a maximum of 159 dB. This adjustment does not 
affect the analog output levels. 

Thus, in order to determine V^g, the nixie tubes or 
the X-Y display are read as follows: 

1. Direct Reading from X-Y Display. 



T = 100 dB 
85 dB 
70 dB 
5 5 dB 
40 dB 



2. Reading from the Nixie Tubes. 



X-Y display 


Nixie tubes 


Remarks 


40 dB 


40 dB 


This reading is present 
when no signal is injected. 


40-80 dB 


40-80 dB 


Actual values on nixies 
are correct 


80-120 dB 


40-80 dB 


Add 40 dB to value shown 
on Nixies to obtain actual 
value 




1" grid 



Calibration: 2" = ^ V 



Reference level 
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The frequencies at which measurements are made are 
displayed both along the x-axis of the X-Y display and on 
the nixie tubes in the following sequence: 



Band 


Frequency 


Band 


Frequency 


Band 


Frequency 


14 


25 Hz 


24 


250 Hz 


34 


2500 Hz 


15 


31.5 


25 


315 


35 


3150 


16 


40 


26 


400 


36 


4000 


17 


50 


27 


500 


37 


5000 


18 


63 


28 


630 


38 


6300 


19 


80 


29 


800 


39 


8000 


20 


100 


30 


1000 


40 


10000 


21 


125 


31 


1250 


41 


12500 


22 


160 


32 


1600 


42 


16000 


23 


200 


33 


2000 


43 


20000 
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APPENDIX C: TERMINOLOGY RELATIONSHIPS 



From Reference 4, the ratio of the EGD probe output 
divided by the corresponding > and then normalized, 
results in 



T' (f) 
g 



I '/An q 

- 8 — ^ 

U °°RMS 



where 



A 

n 



P 

q 



U °°RMS 



spectral probe current 
area of the collector 
number density of charged particles 
charge per particle 

RMS value of the fluctuating component of carrier 
gas velocity. 



Let e = reference voltage of one millivolt. Dividing both 
the numerator and the denominator by e and taking the log- 
arithm of both sides results in 



lo *10 T g (f) = lo «10 ^gio^RMS/e. 

T„ T, 

E h 

where now 

Tg = intensity of the turbulence as measured by the 
EGD probe and expressed in decibels. 

T^ = intensity of the turbulence as measured by the 
hot-wire anemometer and expressed in decibels. 
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